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Abstract Consider the optimal dividend problem for an insurance company whose 
uncontrolled surplus precess evolves as a spectrally negative Levy process. We assume that 
dividends are paid to the shareholders according to admissible strategies whose dividend 
rate is bounded by a constant. The objective is to find a dividend policy so as to maximize 
the expected discounted value of dividends which are paid to the shareholders until the 
company is ruined. Kyprianou, Loeffen and Perez [28] have shown that a refraction 
strategy (also called threshold strategy) forms an optimal strategy under the condition 
that the Levy measure has a completely monotone density. In this paper, we propose an 
alternative approach to this optimal problem. 
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1 Introduction 



The classical optimal dividend problem looks for the strategy that maximizes the expected 
discounted dividend payments until ruin in an insurance portfolio, which has recently re- 
ceived a lot of attention in actuarial mathematics. This optimization problem was first 
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proposed by De Finetti [11] to reflect more realistically the surplus cash flows in an insur- 
ance portfolio, who considered a discrete time random walk with step sizes ± 1 and proved 
that the optimal dividend strategy is a barrier strategy. Since then many researchers have 
tried to address this optimality question under more general and more realistic model as- 
sumptions and until nowadays this turns out to be a rich and challenging field of research 
that needs the combination of tools from analysis, probability and stochastic control. For 
the classical compound Poisson risk model, this problem was solved by Gcrber in [15] via 
a limit of an associated discrete problem. Recently, this optimal dividend problem in the 
classical compound Poisson risk model and also included a general reinsurance strategy as 
a second control possibility was taken up again by Azcue and Muler [7], who used stochas- 
tic control theory and viscosity solutions. For all these cases in general a band strategy 
turns out to be optimal among all admissible strategies. In particular, for exponentially 
distributed claim sizes this optimal strategy simplifies to a barrier strategy. In Albrecher 
and Thonhauser [1] it is shown that the optimality of barrier strategies in the classical 
model with exponential claims still holds if there is a constant force of interest. Avram et 
al. [5] considered the case where the risk process is given by a general spectrally negative 
Levy process and gave a sufficient condition involving the generator of the Levy process 
for optimality of the barrier strategy. Recently, Loeffen [29] showed that barrier strategy 
is optimal among all admissible strategies for general spectrally negative Levy risk pro- 
cesses with completely monotone jump density, and Kyprianou et al. [26] relaxed this 
condition on the jump density to log-convexity. More recent paper Azcue and Muler [8] 
examines the analogous questions in the compound Poisson risk model with investment. 
The corresponding problem in the case of a diffusion risk process was completely solved 
by Shreve et al. [32] and a barrier strategy was identified to be optimal. The special case 
of constant drift and diffusion coefficient was then solved again by slightly different means 
in Jeanblanc-Picque and Shiryaev [21] and Asmussen and Taksar [4]. 

Band strategy (and barrier strategy) often serve as candidates for the optimal strategy 
when the dividend rate is unrestricted. However, the resulting dividend stream is far from 
practical application. In many circumstances this is not desirable. Furthermore, if a band 
strategy is applied, ultimate ruin of the company is certain. Motivated by this fact, other 
dividend strategies such as threshold strategies, linear and nonlinear barrier strategies and 
multi-layer strategies have been studied. Asmussen and Taksar [4] postulated a bounded 
dividend rate and showed that the optimal dividend strategy is a threshold strategy in 
Brownian motion risk models, that is, dividends should be paid out at the maximal 
admissible rate as soon as the surplus exceeds a certain threshold. Some calculations for 
this model can be found in [16]. In the compound Poisson risk model, Gerber and Shiu 
[17] showed that the optimal dividend strategy is a bang bang strategy. In particular, 
for exponentially distributed claim sizes this optimal strategy simplifies to a threshold 
strategy. Motivated by Gcrber and Shiu [17], Fang and Wu [13] study the analogous 
questions in the compound Poisson risk model with constant interest, for the case of an 
exponential claim amount distribution, it is shown that the optimal dividend strategy is 
a threshold strategy. More recently. Fang and Wu [14] examine the same problem for the 
Brownian motion risk model with interest. In a very recent paper, Kyprianou, Loeffen and 
Perez [28] have shown that a refraction strategy (also called threshold strategy) forms an 
optimal strategy under the condition that the Levy measure has a completely monotone 
density. See Albrecher and Thonhauser [2], Avanzi [6] and Schmidli [31] for nice surveys 
on this subject. The purpose of this paper is to re-examine the analogous questions in a 
general spectrally negative Levy process risk model. 

The rest of the paper is organized as follows. In Section 2, we state the problem and 
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recall some preliminaries on spectrally negative Levy processes. In Section 3, we will show 
that the optimal value function of the dividends can be characterized by the Hamilton- 
Jacobi-Bellman (HJB) equation and give a verification result for optimality. In Section 
4 we discuss the threshold strategies. Explicit expressions and the integro-diffcrcntial 
equations for the expected discounted value of dividend payments are obtained, and in 
Section 5 we present the main results, we will show that a threshold strategy forms an 
optimal strategy under the condition that the Levy measure has a completely monotone 
density. 



2 Problem setting 



Suppose that X — {X{t) : i > 0) is a spectrally negative Levy process with probabihties 
{Px : a; G M} such that X(0) = x with probability one, where we write P = Pq. Let 
be the expectation with respect to P^ and write E = Eq. Let {J^t : t > 0} be the natural 
filtration satisfying the usual assumptions. Since the jumps of a spectrally negative Levy 
process are all non-positive, for convenience, we choose the Levy measure to have mass 
only on the positive instead of the negative half line. The Laplace exponent of X is given 

by 

1 

HO) = ae + -a^e^ -J (1 - e-'^ - 9xl{o<.<i}Mdx) (2.1) 

where a e R, cr > and 11 is a measure on (0, oo) satisfying (1 A x'^)Il{dx) < oo and is 
called the Levy measure. The characteristics (a, o"^, 11) are called the Levy triplet of the 
process and completely determines its law. ip is strictly convex on (0, oo) and satisfies 
il!{0+) = 0, -0(00) = 00 and ip'{0+) = EX{1). If cr^ > and 11 = 0, then the process is a 
Brownian motion; When = and Il{dx) < 00, the process is a compound Poisson 
process; When = 0, Il{dx) = 00 and /q°°(1 A x)Il{dx) < 00, the process has an 
infinite number of small jumps but is of finite variation; When o"^ = 0, Il{dx) = 00 
and /q°°(1 a x)Il{dx) = 00, the process has infinitely many jumps and is of unbounded 
variation. In a word, such a Levy process has bounded variation if and only if cr = and 
Jq xll{dx) < 00. In this case the Levy exponent can be re-expressed as 



poo 

m-cP- / {l-e-^nn{dx), 
Jo 



where c = a + xll{dx) is known as the drift coefficient. If > 0,X is said to have a 
Gaussian component. 

In this paper, we shall only consider the case that 11 is absohitcly continuous with 
respect to Lebesgue measure, in which case we shall refer to its density as tt. 

We recall from Kyprianou [24] that for each q > there exits a continuous and 
increasing function W^^^ : R — )■ [0, 00), called the g-scale function defined in such a way 
that W'^'^^x) — for all x < and on [0, 00) its Laplace transform is given by 

/•oo 

/ e-'^W^^\x)dx = , 9 > $(g), (2.2) 

Jo 'H^KP) - Q 
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where = sup{6' > : ip{0) = q} is the right-inverse of ip. We shall write W in place 
of l^^*^) and call this the scale function rather than the 0-scale function. 

The q-scale function takes its name from the identity 



^.(e-^- l(r+ < To") 



iy(«)(a)' 

where 

T+ = inf{t > : X{t) > a}, Tq = inf{t > : X{t) < 0}. 

The following facts about the scale functions are taken from [10, 26]. If X has paths 
of bounded variation then, for all q > 0, W^*^'''*|(o,oo) ^ C^(0,oo) if and only if 11 has no 
atoms. In the case that X has paths of unbounded variation, it is known that, for all 
q > 0, '•'^•'1(0,00) G C^(0, oo). Moreover if cr > then C^(0, oo) may be replaced by 
C^(0, oo). Further, if the Levy measure has a density, then the scale functions are always 
differentiable. In particular, if tt is completely monotone then VF^^)|(o,oo) £ C°°(0, oo). It 
is well known that W^^\0+) = 1/c when X has paths of bounded variation. Otherwise 
iy(^)(0+) = for the case of unbounded variation. In all cases, if EX[1) > 0, then 
W{oo) = 1/EX{1). If g > 0, then W^i\x) - e*('')^/V^'($(g)) as a; ^ oo. 

Spectrally negative Levy processes have been considered recently in [12, 18, 20, 23, 
25, 34], among others, in the context of insurance risk models. It is assumed that, in 
the absence of dividends, the surplus of a company at time t is X{t). We assume now 
that the company pay dividends to its shareholders according to some strategy. Let 
^ = {Lf : i > 0} be a dividend strategy consisting of a left-continuous non-negative 
non-dccrcasing process adapted to the filtration {J-'t : t > 0} of X. Lf represents the 
cumulative dividends paid out up to time t under the control ^ by the insurance company 
whose risk process is modelled by X. We define the controlled risk process — {U^i^) : 
i > 0} by U^{t) = X{t) - Lf . Let T = inf{i > : U^{t) < 0} be the ruin time, and define 
the value function of a dividend strategy ^ by 

V^{x) = Q\-i'dLiy 

where g > is the discounted rate. 

A dividend strategy is called admissible if Lj_,_ — I/f < U^{t) V for t < T, in other 
words the lump sum dividend payment is smaller than the size of the available capitals. 
Let S be the set of all admissible dividend policies. The control problem consists of solving 
the following stochastic control problem: 

K(x) = sup V'^(x), 

and, if it exists, to find a strategy G S such that V^*{x) = for all x >0. 

In this paper, wc assume that the admissible dividend rate is r(t) at time t which is 
bounded by a constant a. In the sequels, we assume that < a < a + xn{x)dx if X 
has paths of bounded variation. Under this additional constraint, we will show that if 
the Levy measure 11 has a completely monotone density and that 5 > 0, then the optimal 
dividend strategy is formed by a threshold strategy. 
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3 The HJB equation and verification of optimality 



Let 

V{x) = sup Vr{x), 
r(.) 

where the supremum is taken over all control process r(t) which are admissible according 
to the constraints. Vr is the value function when the admissible dividend rate is r{t) 
at time t and V is called the optimal value function. Suppose V is twice continuously 
differential on (0, oo) when the process X is of unbounded variation and is continuously 
differential on (0, oo) when the process X is of bounded variation. Standard Markovian 
arguments yield that V satisfies the following Hamilton- Jacobi-Bellman (HJB) equation 
(the proof is similar to the one in Azcue and Muler [8]): 

max [1 - V'{x)]r + TV{x, h) - 5V{x, 6) = 0, x > 0, (3.1) 

0<r<a 

where 

VV{x,h) = \a^V'\xM + aV\x,h) 

+ j^[V{x -y,b)- V{x, b) + V'{x, b)yl^o<y<iMy)dy. 

From above expression, as in Gerber and Shiu [17], we see that at time t G (0,T), the 
optimal dividend rate is 

r = if V'{U^{t-)) > 1, 

r = a if V'{U^{t-)) < 1. 

If V'{U^(t—)) = 1, the dividend rate r can be any value between and a. In particular, 
the optimal dividend rate at time is 

r = if V'{x) > 1, 

r = a if V'{x) < 1. 

Thus, the company should either pay nothing, or the maximum possible. This called a 
bang bang strategy. 

Now we show that a strategy ^ with z/(x) = V^{x) is smooth enough and satisfying 
the HJB equation (13.11) is indeed an optimal strategy. 

Consider any other dividend strategy, with dividend rate r(t) and surplus X(t) at time 
t. We claim that ^ 

E (^j^ e-^'r{t)dt\X{0) = x^ < u{x). (3.3) 

From this, it follows that z/(x) = V{x), and hence the given strategy ^ is optimal. 
To prove inequality (13.31) . we consider the martingale 

e-^'u{X{t)) - e-^' ((F - 6)u{X{s) - r(s)z/'(X(s)) ds, 

5 



which can be shown by Ito's formula for semimartingale. From optimal sampling theorem, 
we have 

E (^e-'^(*^^V(X(t AT) - e'^' (^{T - 5)u{X{s)) - r(s)z/'(X(s)) cis|X(0) = = u{x), 



which implies 



E (^j^ e~^' ((r - 5)u{X{s)) - r(s)z/'(X(s)) rfs|X(0) = x^ < u{x) 



smce 



E (e-'^(*^'^)z/(X(t A T))|X(0) = > 0. 
Because the function z/(x) satisfies the HJB equation (13. ip . we have 

r(s) + (r - 6)u{X{s)) - r{sy{X{s)) < 0. 

Thus 



tAT 

5s 



E / e-''r{s)ds\X{0) = x 



tAT 



- S)uiX{s)) - r{s)i^\X{s))ds\X{0) =x] < u{x). 

Letting t — > oo yields (13.31) . 



4 Threshold dividend strategies 



In this section, we study the threshold strategy. We assume that the company pays 
dividends according to the following strategy governed by parameters 6 > and a > 0. 
Whenever the modified surplus is below the threshold level b, no dividends are paid. 
However, when the surplus is above this threshold level, dividends are paid continuously 
at a constant rate a that does not exceed the premium rate c. Note that if a = c, we 
have a barrier strategy again. We define the modified risk process Uh = {Ub{t) : t > 0} by 
Ub{t) = X{t) — Dh{t), where -Db(t) = a l{Ub{t) > b)dt. The existence of such a process 
and some conclusions on fluctuation identities can be found in Kyprianou and Loeffen 
[27]. Let Dfy denote the present value of all dividends until time of ruin T, 

Db = a [ e-^'l{Ub{t) > b)dt 
Jo 

where T = inf{t > : Ub{t) < 0} with T = oo if Ub{t) > for all t > 0. Here 5 > is the 
discount factor. Denote by V{x, b) the expected discounted value of dividend payments, 
that is, 

V{x,b) = E{Db\Ub{0) = x). 
Clearly, < V{x, 6) < f and lim^^oo V{x, b) = f . 
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Define the first passage times, with the convention inf = oo, 

T+ = inf{t > : Ub{t) > 6}, T" = inf{t > : Ub{t) < b}. 

Let Y = {Y{t) := X{t) - at}t>o. For each 6 > 0, W^^^ and Z^^^ are the 5-scale 

functions associated with X and wi^^ and z!f^ are the (5-scale functions associated with 
Y. Further, \1/ is defined as the right inverse of the Laplace exponent of Y so that 

^{6) = sup{e > : ipie) -ae = 6}. 

Theorem 4.1. Assume W^^^ is continuously differentiable on (0,oo). 
(1) For < X < b, we have 

n-,'')-V(b,i)'^. (4.1) 



(2) For X > b, we have 

HS) 



V{x,b) = -a j;-' W^'\y)dy + ^W^'\x - b) 



)(b) 

where 



+^V{b,b) (W^^y{x - 6) - '^{6)Wy\x - 6)) (4.2) 

r /-oo«i'^(^ -b,y- bMv - z)w('H^)dydz. 



u 



x,y) = Wi'\x)e-''^'^y -Wi'\x-y). 



Proof (1). For < x < 6, using the strong Markov property of U at T^"*", we have 

V{x,b) = y(6,6)E,(e-<l(T+ < T)), 

and (14. ip follows since 

E^{e-'''^' 1(T+ < T)) 



W(^){b)' 

(2). For X > b, using the strong Markov property of Ub at T^~, we have 

V{x,b) = ^P,{T,~ = oo) + aE, (^f^^~ e-'^du,T,- < oo 
+aE^ J^- l{Ub{u) > b)e~^''du,T- < oo) 

= fP.(T,- = oo) + fE.((l-e-^^^"),T,-<oo 
+E, (e-'''^~V{Ub{T,-),b),T,- < oo 

= ^ - jE^ {e-'"^^' < oo) 

+E^ (e-'^^V{Ub{T^),b).Tb < oo 



(4.3) 
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Note that 



E^{e-''">'ViUtiT-),b),T-<oo 



V{h,h)E^(e-'^^ MT-) =h) 

poo pb 

+ / ui^\x - b,y - b)n{y - z)V{z,b)dydz 

2 

yV(6,6) (w,^^'^\x - 6) - ^i5)Wi^\x - hi 

POO pb 

+ / ui^\x - b,y - b)n{y - z)V{z,b)dydz. 

Jb J — oo 



In particular, 



where 



It follows that 



V{x,b) 



y (W.^^'^'ix - 6) - )(x - b) 

r>00 

/ u^^^ {x-b,y - b)'K{y - z)dydz 
Zi'\x-b)-^^Wi'\x-b), 



Zi'\x) = l + 6 / W^'\y)dy. 
Jo 



2 



+ 



oo pb 



uf^ {x-b,y- b)7r{y - z)V{z, b)dydz. 



Putting (14. ip into above expression leads to f l4.2p . 

The following result agrees with the result of Kyprianou and Loeffen [27, (10.25]. 

Corollary 4.1. Suppose X has paths of hounded variation and let < a < c, where 
c = a + xll{dx). 

(1) For < X < b, we have 

V{x,b) 

(2) For X > b, we have 

, , ,,,, , W^^^ ^aCW^^\x-y)W^^^'{y)dy , ^ 



X) 



^{5)e^i5)b e-'^(^>W(^y{z)dz 



(4.4) 



Proof (1). It follows from Kyprianou and Loeffen [27, (4.10] that, for x > b, 

POO pb 

/ / u'^^\x-b,y-b)n{y-z)W^^\z)dydz 

J b J — oo 

)W^^^\z)dz 

POO 



Substituting this into (14 ■2p and letting x — )■ 6 we find that 

W^^\b) 



V{b,b) 



This, together with (gj]) and (i2]), leads to (143D and f H3|) . 

The next result was obtained by Gerber and Shiu [17] for the compound Poisson 
model: 

Theorem 4.2. Suppose that X has no Gaussian component. Then, as a function of x, 
V{x, b) satisfies the following integro- differential equations: 

aV'{x,b) + J^°^[Vix-y,b)-Vix,b) + V'ix,b)yl^o<y<iMy)dy ^^^^ 
= 5V{x, b), < X <b, 



{a-a)V'{x,b) + j;;^[V{x -y,b) -V{x,b) + V'{x,b)yl^o<y<i)]7riy)dy ^^^.^ 
= 6V{x, b) — a, X > b, 

with the continuity condition V{b—,b) = V{b+,b) = V{b,b). Moreover, if X has paths of 
bounded variation then 

cV'ib-, b) = {c- a)V'{b+, b) + a; 
// X has paths of unbounded variation then 

V'{b-,b) = V'{b+,b). 

where c = a + xT^{x)dx. 



Proof Equations (14.61) and (14. 7p can be proved by Ito's formula. It follows from (14.11) 
that V(6— ,6) = V^(6, 6), since W^^"* is continuous at b. From fl4.2p we see that V{x.,b) is 
differential on [6, oo), where we at x = 6 mean the right-hand derivative. Consequently, 
V{b+,b) = V{b,b). This proves the continuity of V at b. If X has paths of bounded 
variation, then the derivative, V'{x,b), is not necessarily continuous at x = 6. In fact, it 
follows from equations (14.61) and (14. 7p that 

cV'{b-, b) = {c- a)V'{b+, b) + a. 
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If X has paths of unbounded variation, let n„ be measures on {1/n, oo): 

n„(rfx) = Il{dx)l(i/n,oo), n>l. 

Then xlln{dx) < oo. This is to say that a process X„ with the Levy measure n„ has 
paths of bounded variation. From above we get 

c„K(6-, b) = (c„ - a)V;{b+, h) + a, 

where Cn = a + XTTn{x)dx. Letting n — t- oo yields 

V'{b-,b) = V'{b+,b). 

This ends the proof of Theorem 4.2. 

The next result was obtained by Wan [33] for the compound Poisson model perturbed 
by diffusion: 

Theorem 4.3. Suppose that X has a Gaussian component a > 0. Then, as a function 
of X, V{x, b) satisfies the following integro- differential equations: 

id b) + aV'ix, b) + j^[V{x -y,b)- V{x, b) + V\x, 6)i/l(o<,<i)]7r(|/)rfy 
= 5V{x, b), < X <b, 

(4.8) 

'^a'V'ix, b) + {a- a)V'ix, b) + j^[V{x -y,b)- b) + V'{x, b)yl^o<y<iMy)dy 
= 6V{x, b) — a, X > b, 

(4.9) 

with the boundary conditions V{0,b) = 0. Moreover 

V{b-,b) = Vib+,b) = Vib,b), V'ib-,b) = V'{b+,b). 

Proof Equations (14.81) and (14. 9 p can be proved by Ito's formula. If Ub{0) = 0, because 
cr > 0, ruin is immediate and no dividend is paid, so we have V^(0, b) = 0. It follows from 
(14. ip that V{b—,b) = V{b,b), since W^^^ is continuous at b. From (14. 3 p we have 

V{b+, b)<j(^l- E,e-'^^^ + V{b, b)E, [e-'^^'^ = V{b, b), 

and thus V{b+, b) = V{b, b). This proves the continuity of V at b. 

In the following, we prove that {aB{t);t > 0} can be approximated by the process 
{cet — eN^(t); t > 0}, where N^(t) is a Poisson process with parameter > 0, and 
Ce > is a constant. Now, we choose e, A^, and such that Var[eN^(t)] = aH and 
E[cet — eN^{t)] = 0. These two conditions yield A^ = cr^/e^ and = a'^/e. It is easy 
to prove that, when e — ?► 0^, _E'[e^(^=*^^^-(*))] " This shows that the process 

{cet — eNe(t); t > 0} converges weakly to the process {crBit); t > 0}. It follows that the 
Levy process X with Levy triplet (a, a, 11) can be approximated by the Levy process 
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with Levy triplet (a + 0, Ilg), where Ilg = 11 + l(a;>£). Therefore, by Theorem 4.2, for 
example, in the case of bounded variation, we have 

(c + c,)V^{b-, b) = {c + c,- a)V^{b+, h) + a. 

Letting £ — >■ and noting that c, a, V^'(&— , h) and V^(&+, &) are bounded, and lim^-^o K = 
V , yields 

V^'(6-,6) = V'{h+,h). 
This ends the proof of Theorem 4.3. 

5 Optimal dividend strategies 

In some situations the optimal dividend strategy is a threshold strategy. It is easy to 
see that if V'{x^ 0) < 1 for a; > 0, then the threshold strategy with 6* = is optimal, if 
V'{x,h*) > 1 for x < 6* and V'{x,h*) < 1 for a; > 6*, then the threshold strategy with 
6* > is optimal. The optimal threshold h* can be obtained V'(h*,h*) — 1. In fact, for 
the case cr > 0, it is clearly, since V'{x, h) is a continuous functions of x on (0, oo); For the 
case o" = 0, using the the same argument as in Gerber and Shiu [17], the result follows. 
From those facts one sees that if ^(a;, 6*) is a continuously differentiable concave function 
on (0, oo), then the optimal dividend strategy is a threshold strategy. 

We now review definitions and some properties of logconvex functions and completely 
monotone functions. We refer the readers to [3, 9] for more details. 

A function / defined on an convex subset of a real vector space and taking positive 
values is said to be logarithmically convex if log(/(x)) is a convex function of x. It is easy 
to see that a logarithmically convex function is a convex function, but the converse is not 
always true. For example f{x) = is a convex function, but log(/(x)) = 2 log \x\ is not 
a convex function and thus f{x) — x^ is not logarithmically convex. 

Recall that a / G C"^(0,oo) with / > is completely monotone if its derivatives 
alternate in sign, i.e. (— l)"/^ > for all n e N. 

Note that the class of logconvex functions contains the class of completely monotone 
functions. In fact, any completely monotone function is both nonincreasing and logconvex. 

Some distributions with completely monotone density functions are (see [9, 29]): 

• WeibuU distribution with density: j{x) = crx'^~^e~'^^'^ , x > 0, with c > and 
< r < 1. 

• Pareto distribution with density: f{x) — a{l + x)'"''^, x > 0, with a > 0. 

• Mixture of exponential densities: f{x) = ^^Li ^iAe^^'"*, x > 0, with Ai > 0, (3i > 
for i = 1, 2 • • • ,n, and Yli=i — 1- 
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• Gamma distribution with density: f{x) = ^-fj§^0--, a; > 0, with /3>0,0<c<l. 

The following are several important examples of spectrally negative Levy processes 
with completely monotone densities and that satisfy Tc{x)dx = oo (cf. [19, 29]): 

• a-stable process with Levy density: 7c{x) = Ax~^~", x > with A > and 
ae (0,1)U(1,2); 

• One-sided tempered stale process (particular cases include gamma process (a = 0) 
and inverse Gaussian process (a = ^)) with Levy density: n^x) = Ax~^~"e~'^^, x > 
with /3, A > and -1 < a < 2; 

• The associated parent process with Levy density: 7r(x) = Xix~^~°'e~^^+X2X~'^~°'e~^^ , 
X > with Ai, A2 > and — 1 < a < 1. 

More examples can be found in recent paper of Jeannin and Pistorius [22, Example 
2.4]. 

The following result can be found in Loeffen and Renaud [30]: 

Lemma 5.1. Suppose the tail of the Levy measure is log-convex, then, for all 6>0, W^^^ 
has a log- convex first derivative. 

Theorem 5.1. Suppose the tail of the Levy measure is log-convex, then, for all 6 > 0, 
V{x,b*) is a concave function on (0,6*), where b* is the solution of V'{b* ,b*) = 1. 



Proof Differentiate (14. ip with respect to x, and then set x = b* yields 

V'(x, b*) = ,,,, ' , 0<x<b*. (5.1) 

^ ^ w^^y{b*) ~ ~ ^ ^ 

Since b* is the value of b that maximizes V{x, b), i.e. b* is the value where W^^^i^b) attains 

its global minimum, and thus W^^^ (x) is decreasing on (0, b*). It follows that V'{x, b*) is 
decreasing on (0,6*), which implies V{x,b*) is a concave function on (0,6*). 

Theorem 5.2. Suppose that the Levy density it is a completely monotone function on 
(0,00) and that 6 > 0. Then V{x,b*) is a concave function on (6*, 00), where 6* is the 
solution ofV'{b*,b*) = 1. 



Proof We prove the theorem in three cases. 

Case 1: 11(0, 00) < 00 and a = 0. Using (14.51) and repeating the proof of Lemma 8 in 
Kyprianou, Loeffen and Perez [28] we get V"{x,b*) <0,x > 6*. 

Case 2: 11(0, 00) = 00 and cr = 0. By Bernstein's theorem we can express the function 
TT in the form: 

- 00 



7r(x) = / e """'//((iM) 
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where is a measure on (0, oo). It is well known that tt can be approximated arbitrarily 
closely by a hyper-exponential density (by approximating the measure /j, by sums of point 
masses). It follows from Jeannin and Pistorius [22] that X can be approximated by a 
sequence of approximating processes (X'^"^)„>i. The approximating density equal to 

i 

where {ui)i = {u^P^)i is the finite partitions of (0, oo), and (Aj)j = (A|"''*)i is finite sets 
of positive weights. For a given n, the partitions {ui)i and the weights (Ai)j are satisfy 
certain conditions. The approximating process X*^"^ constructed in this way can be shown 
to converge weakly to X. For more details, see Jeannin and Pistorius [22, §3.2]. For each 
n > 1, we find that 7r„(x)dx < oo and 7r„ is complete monotone on (0, oo). Let Vn be 
the expected discounted value of dividend payments corresponding to the Levy process 
{X^"'^)n>i. From the Case 1, we know that V^(x, 6*) is a concave function on (6*,oo), 
where 6* is the solution of V^{b'^,bD — 1. The result follows, since lim„^oo V„ = V, 
lim„^oo ^n — ^' ^iid the limit of a pointwise convergent sequence of concave functions is 
concave. 

Case 3: cr 7^ 0. We consider the following approximating process (X('*))„>i (cf. Loeffen 
and Renaud [30]): The Levy triplet is (an,0,n„), where a„ = a + |(T^ne^"'(n + 1) and, 
n„ is defined by 

n„(x, 00) = n(x, 00) + -cr^n^e""^. 

For all ^ > 0, lim„_^oo V'nC^) = '0(^)) thus, by the continuity theorem for Laplace 
transforms, hm„_^oo — W'^^^ {x) for all x > 0. Let Vn be the expected discounted 

value of dividend payments corresponding to the Levy process X("). It is easy to see that 
n„ has a completely monotone density and thus from the Case 1 and Case 2, we know 
that 6* ) is a concave function on (6* , 00), where 6* is the solution of V^^(&^, — 1. 
Because lim„^oo Vn — V, hnin-^-oo V^ — V, the result follows. This completes the proof of 
Theorem 5.2. 

Combining Theorems 5.1 and 5.2 we get the main result of this paper which established 
by Kyprianou, Loeffen and Perez [28] by using an alternative argument: 

Theorem 5.3. Suppose that the Levy measure U has a completely monotone density and 

that S > 0. Then V{x,b*) is a concave function on (0, 00), where b* is the solution 
of V'{b* ,b*) = 1. Consequently, the threshold strategy with threshold b* is the optimal 
dividend strategy. 

Remark 5.1. If there were no restrictions on the admissible dividend rate r(t), or, a = c 
in the case of X has paths of bounded variation, then the threshold becomes a barrier. 
Suppose that the Levy measure 11 has a completely monotone density and that 6 > 0, 
then barrier strategy at b* is an optimal strategy. This result can be found in Loeffen 
[29] . Kyprianou et al. [26] strengthen the result of Loeffen [29] estabhshed a larger class 
of Levy processes for which the barrier strategy is optimal among all admissible ones. An 
alternative proof is given in Yin and Wang [35] . 

Remcirk 5.2. We can conclude that the threshold strategy is the optimal dividend strat- 
egy for the compound Poisson risk model or the compound Poisson risk model perturbed 
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by Brownian motion where the claims have a distribution with a completely monotone 
probability density function, which extended the known result (c.f. Gerber and Shiu [17] 
where only consider the exponential claim density under the Cramer- Lundberg setting). 
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